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SPECTRAL FUNCTIONS RELATED TO SOME FRACTIONAL 
STOCHASTIC DIFFERENTIAL EQUATIONS 

MIRKO D’OVIDIO, ENZO ORSINGHER, AND LUDMILA SAKHNO 


Abstract. In this paper we consider fractional higher-order stochastic differential equa¬ 
tions of the form 

+ f > 0, ^ > 0, ^ > 0, a € (0,1)UN 

where £{t) is a Gaussian white noise. We derive stochastic processes satisfying the 
above equations of which we obtain explicitly the covariance functions and the spectral 
functions. 


1. Introduction 


In this paper we consider fractional stochastic ordinary differential equations of different 
form where the stochastic component is represented by a Gaussian white noise. Most of 
the fractional equations considered here are related to the higher-order heat equations and 
thus are connected with pseudoprecesses. 

The first part of the paper considers the following stochastic differential equation 



X{t) 


£{t), 


/3 >0, 0 < a < 1, 0 


( 1 . 1 ) 


where represents the upper-Weyl fractional derivative. We obtain a representation 

of the solution to (1.1) in the form 

1 POO POO 

X{t) = —— / dz j dss^~^e~''^ha{z,s)£{t +z) (1.2) 

1 (P) Jo Jo 

where ha{z, s), z, s > 0, is the density function of a positively skewed stable process Ha{t), 
t > 0 of order a. G (0,1), that is with Laplace transform 

e~^^ha{z,s)dz = , ^> 0 . 

For (1.2), we obtain the spectral function 



fi-r) 


(1.3) 


and the related covariance function. 

The second type of stochastic differential equations we consider has the form 




Q2n 


P 


X{t) = £{t), P > 0, fj. > 0, n> 1, 


(1.4) 
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where £{t) is a Gaussian white noise. The representation of the solution to (1.4) is 

1 f'oo r+oc 

X{t) = —^ / dxu 2 n{x,w)£{t + x) (1.5) 

1 [P) Jo J-oo 

where U 2 n{x,w), x G M, tt; > 0 is the fundamental solution to 2 n-th order heat equation 

-{x,w) = {-lT^^^{x,w) ( 1 . 6 ) 

The autocovariance function of the process (1.5) can be written as 

^2 t-oo 2 

EX{t)X{t + h) =—— dww'^^~^e~^'^ U 2 n{h,w) = ^Eu 2 n{h,W 2 is) (1.7) 

1 y^P) Jo P ^ 

where W 2 j 3 is a gamma r.v. with parameters /r and 2/3. The spectral function /(r) 
associated with (1.7) has the fine form 

/M = +‘(-2n)28 (1'8) 

For n = 1, (1.6) is the classical heat equation, U 2 {x,w) = -^ 71 = and, from (1.7) we obtain 
an explicit form of the covariance function in terms of the modified Bessel functions. 

In connection with the equations of the form (1.6) the so-called pseudoprocesses, first 

introduced at the beginning of the Sixties ([ 6 ]) have been constructed. The solutions to 
( 1 . 6 ) are sign-varying and their structure has been explored by means of the steepest 
descent method ([9; 1]) and their representation has been recently given by [10]. 

For the fractional odd-order stochastic differential equation 

(^/z + (-l)"^^j X{t)=£{t), n = l,2,... (1.9) 


the solution has the structure 


X{t) 


1 

m 



dx U 2 n+l{x, w)£{t + x) 


where U 2 n+iix,w), x G M, u) > 0 is the fundamental solution to 


( 1 . 10 ) 


The solutions U 2 n-i-i and U 2 n are substantially different in their behaviour and structure 
as shown in [10] and [7]. 

A special attention has been devoted to the case n = 1 for which (1.10) takes the 
interesting form 




1 



dx 




£{t x) 


( 1 . 12 ) 


where Ai{-) is the first-type Airy function . The process X 3 can also be represented as 

X^{t) = ^E£{t + Y^{Wp)) (1.13) 

where the mean value must be meant w.r.t. ^^(IF^g) and I 3 is the pseudoprocess related 
to equation 

du _ d^u 
dt dx^ 


(1.14) 
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and Wjs is a Gamma-distributed r.v. with parameters 13, ji, independent from Y^. The 
autocovariance function of has the following form 


¥.X^{t)X^{t + h) = 


1 ^7 h y 


(1.15) 


where VF 2/3 is the sum of two independent r.v.’s Wj^. For the solution to the general 
odd-order stochastic equation we obtain the covariance function 

2 

EX{t)X{t + h) = ^E [U2n+i{h, W2p)] (1.16) 

Of course, the Fourier transform of (1-16) becomes 

Stochastic fractional differential equations similar to those dealt with here have been 
analysed in [2], [3] and [5]. In our paper we consider equations where different operators 
are involved. 


2. A STOCHASTIC EQUATION INVOLVING FRACTIONAL POWERS OF FRACTIONAL 

OPERATORS 

In this section we consider the following generalization of the Gay and Heyde equation 
(see [3]) 

/i -I- ^ ) X{t) = S{t), /3>0, 0<q:<1, /j. > 0 


d{-tY 

where T(t), t > 0, is a Gaussian white noise with 


E£{t)£{s) = 


a‘‘, t = s 


0 , t Y s 

The fractional derivative appearing in (2.1) must be meant, for 0 < a < 1, as 




d{-ty 


Jit) = - 


1 


d 


_^ fj) 

F(1 — a) dt {s — ty 


-ds 


a 


F(1 - a) do 


f 


fit) - fit + w) 


w 


a+1 


-dw. 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 


For information on fractional derivatives of this form, called also Marchaud derivatives, 
consult [11, pag. 111]. For A > 0, we introduce the Laplace transform 


C 


d^f 

[d(=t)»\ 

which can be immediately obtained by considering that 




e-w^^/(t)dt = A“£[/](A) 


£ 


d'^f 


(A) = 


a 


(£|/l(A)-e-“V[/|(A)) 


dw 


w 


a+l 


(2.5) 


( 2 . 6 ) 


F(1 - a) Jo 

for a function / such that e^^f{t) G L^([0,oo)). 

Theorem 2 . 1 . The representation of a solution to the equation (2.1) can be written as 

1 


^00 poo 

X{t) = —J— / dz j dss^~^e~‘'^haiz,s)£{t +z) 
1 yp) Jo Jo 


(2.7) 
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Proof. The solution to the equation (2.1) can be obtained as follows 


Xit) = 


roo 


d{-t) 

mJo 


+ h £{t) 


— I —Sfl—S 


r(/3) JO 

1 


ds 




e } ds. 


r(/3) Jo 

Now, for the stable subordinator t > 0, we have that 

roo ^ 

= / dzha{z,s)e^^£{t) 

Jo 

poo 

= / dz ha{z, s) £{t + z) 


( 2 . 8 ) 


(2.9) 


where ha{z, s) is the probability law of H°^{s), s > 0. In the last step of (2.9) we used the 
translation property 


e^dt£{t) = £{t + z). 


This is because 


z — 

e dt 


OO h ih 

/yhj 

k=0 


In view of the Taylor expansion 

OO / \k 

f(x) = 


k=0 


with xq = t and x = t + z we have that 


e^dt = f>{t + z) 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


which holds for a bounded and continuous function cj) : [0, oo) i—)• [0, oo). Since we can find 
a sequence of r.v.’s {ojjjgN and an orthonormal set, say for which (2.13) holds 

true V j and such that 


lim E 

N—^OO 


N 






i=i 


= 0 , 


we can write (2.10). Therefore, 

1 poo poo 

X{t) = —— / dz ds haiz, s) £{t + z) (2-14) 

-1 (P j Jo Jo 

is the formal solution to the fractional eqnation (2.1) with representation, in mean sqnare 
sense, given by 


iSN 


□ 


Remark 2.1. For the case a = 1, haiz, s) = 6{z — s) where 6 is the Dirac delta function 
and from (2.10) we infer that 


e ^^s^ ^ £{t + s) ds 


(2.16) 
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solves the fractional equation 


X{t) =8{t). 


(2.17) 


Consult on this point [5]. 

A direct proof is also possible because from (2.8) we have that 


Xit) = 


/-I (i) ds 


r(/3) Jo 

1 

^ Jq + (2.18) 

In the last step we applied (2.10). 

Remark 2.2. For a = 1 and (3 = 1, we observe that (2.7) becomes the Ornstein-Uhlenbeck 
process. 

Our next step is the evaluation of the Fourier transform of the covariance function of 
the solution to the differential equation (2.1). Let 

r‘+0O 


where 


/ + 00 

e^^^Covx{h)dh 

■OO 

Covx{h) =E[X{t + h)X{t)] 


with EX(t) = 0. 

Theorem 2.2. The spectral density of (2.7) is 


Ht) = 


a 


{pf + 2|r|"//cos ^ + |r|2“) 


, r G M, 0 < a < 1, /3 > 0. (2.19) 


Proof. The Fourier transform of the covariance function of lag h = t 2 — ti of (2.7) is given 
by 

roo poo 

/ / e"^^^-^^^EX{ti)X{t2)dtidt2 

Jo Jo 

1 POO POO POO POO POO POO 

^^2 / dzi / dsi / ds 2 / dz 2 S^~^S 2 ~^ 

[PI Jo Jo Jo Jo Jo Jo 


F2(/3) 

X {zi, Si) ha {Z2, S 2 ) E£{ti + zi)£{t2 + Z 2 ) 


where 


Thus, 

r*oo poo 


E£{ti + zi)£{t2 + Z2) = 


(T^, h = Zl - Z2 
0, h ^ Zl - Z2 


( 2 . 20 ) 


POO POO 

/ / e*^(*2-L)EX(ti)X(t2)dtidt2 = 

Jo Jo 


fj" 

Jo 


OO poo poo poo 


dzi / dsi / ds2 / dz2 sf ^ ^ 

Jo Jo Jo 

X (zi,Si) ha (Z2,S2) 

By considering the characteristic function of a positively-skewed stable process with law 
ha, we have that 


e^'^^^ha{zi,si) dzi = e 


/ • NCK I ifv —^-^T sqriT 

^ p-(-*T) SI _ g-si|r|“e ^ 


and 


e /la {Z 2 , S 2 ) dz 2 = e 


^ ^-{irrs2 ^ g-S2|rre*^=«""_ 


( 2 . 21 ) 

( 2 . 22 ) 
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Thus, we obtain that 

r*oo /•oo 


roo roo 

Jo Jo 


^2 /•oo /•oo 

^ ’ /5-1 „/3-1„-(si+S2)m p-(*'r)“s2-(-iT)“si 


p(/3)7„ 


a 


(/X + |r|“ e-^ ^snry y ^ sgnry 


a 


(/x^ + 2|r|“/xcos ^ + |r|2") 


/?• 


□ 


Remark 2.3. In the special case a = 1 the result above simplifies and yields 

,2 


/('t) = 


a 


(/x2 + ry 


(2.23) 


We note that for fi = 1, (2.23) becomes the spectral function of the Ornstein-Uhlenbeck 
process. Processes with the spectral function f are dealt with, for example, in [2] where also 
space-time random fields governed by stochastic equations are considered. The covariance 
function is given by 


Covx{h) = 7 ^ [ e *^^/(r)dT 

\m Jo 

a/ 

Jo Jk 


/ e—" 

271 JjR 
^2 poo 


m 

^2 poo 


) dr 


) dz 


a 


r(/?) 


2r(/3)r( 


poo 

/ z^-^ 

Jo 

of 

g) Jo 


2 6 42 

e T^dz 
y/Airz 




cr 


2 /ILI\/3-| 


\h\ 


r(/?)r(i) V 2 /uy 1 ^’ - 

where Ky is the modified Bessel function with intergal representation given by 

j exp {—/3x^ — dx = - Ke , p,a,fi,v>t) (2.24) 


Ky{x) « - for X —)■ 0"' 


(see for example [4], formula 3.478). We observe that Ky = K^y and Ki{x) = ^■ 

(2.25) 

(2.26) 
(2.27) 


Moreover, 

([8, pag. 136]) and 

Thus, we get that 


Kyix) K. . —e ^ for x —)> oo. 
2x 


Covx{h) ~ for /i —)• 0^ 
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and ^ 

Covx{h) ^ /i —5-oo. (2.28) 

We study the covariance of (1.2). Recall that, a stable process S of order a with density 
g is characterized by 

g{C, t) = a G (0, 2]. 

Consider two independent stable processes Si{w), S 2 {w), w > 0, with = 1 and (t| = 
2/xcos^. Let gi{x,w), x G M, rc > 0 and g 2 {x,w), x G M, tc > 0 be the corresponding 
density laws. Then, the following result holds true. 

Theorem 2.3. The covariance function of (1.2) is 

^2 roo r-\-oo 

Covx{h) = :r:^ gi{h-z,w)g 2 iz,w)dzdw (2.29) 

-L [P) Jo J-oo 

or 

2 

Covx{h) = -^Egs^+S 2 ih,Wis) (2.30) 

and is a gamma r.v. with parameters g‘^.,13. 

Proof. Notice that 

^/3-lg-«;(M2+2|r|“MCos 

where 

^- 2 ^^cosr^\T\-w ^ ^^iTS 2 {w) ^ and = gi{T,w). 



Thus, 


a 


h 

from which, we immediately get that 


/(r) = —E[5i(r,W^)52(r,VL/3)] 


a 


CovX{h) = 7 ^E 


P 


r-+00 


gi{h - z,Wp)g 2 {z,W^)dz 


a 


r(/3) 


/•OO 

^0 


r-\-oo 




gi{h — z, w)g 2 {z, w)dz dw 


□ 


3. Fractional powers of higher-order operators 
We focus our attention on the following equation 



X{t) = S{t), 


t>0, ^>0, /3>0 


(3.1) 


that is, to the equation (1.4) for n = 1. 

Theorem 3 . 1 . The representation of a solution to the equation (3.1) can be written as 

roo ^+c>o 

X{t) = / w^~^e~^'^ / U 2 {x,w)£{t + x) dx dw, /3 > 0, g>0. (3.2) 

r (p) Jo J-oo 

Moreover, the spectral function of (3.1) reads 


f{T) 


(p r2)2/3 


(3.3) 
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and the corresponding covariance function has the form 


a 


Covx{h) 




e 4^2,3 




2cj" 


h\_\ 

r(2/3) \ 2 ^) 


2/3 


K2p{\h\^) 


(3.4) 


where W 2/3 is a gamma r.v. with parameters g,, 2/3. 
Proof. We can formally write 


e = 


f 


e 4 ™ 

e dt _ dx 


2y/'KW 


(3.5) 


so that from (3.1) we have that 

X{t) = ^ 


r{P)J„ 

/*oo 


e-^^wf^- 


r(/3) 


/ e 

Jo 


-M«',„/3- 


i-c 

Uw / 

J —C 

/ C 

-c 


TTU’ 

2 


e^®f(t) dx 


TTIC 


£{t + x) dx. 


(3.6) 


By observing that 


we can write 




(T^, Xl — X2 = h 


^2 poo POO P 

EX{t)X{t + h) =Y2-^ J e~^'^^w^~^dwi J e~^'"^w^~^dw2 J 


POO POO 

/ e-^^^w^-^dwx / e-^^^w{-^dw2—r 

Jo Jo 2 / 


^ /p 4w 


e e 


jh-xiy 

4wo. 


-00 2 ^/ 7 raI 2^/7™^ 


-dxi 


a 


g 4(tD2+m2) 


r2 (/3) 

2 


7r(t(;i + 7x2) 








g 4(vri+W2) 


e 

2 VJJW 


a 


r (2/3) do 20™ 
2^2 ( h 


rx 


2/3-ig-M«^dR; 


r(2/3) V20I 


(—)■ 

V2VM/ 


K2l3{hy/Ji) 


We notice that 


u 


Covx{h) = ^P{B{W 2 g) G dh)/dh 


P 


where B{W 2 g) is a Brownian motion with random time W 2 / 3 . Thus, we obtain that 


/ OO 

P-’^Covxih) 

■00 


dh = 


/t2 poo 

Jo 


a 


r(2/3) 


(^ + r2)2/3' 


□ 
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An alternative representation of the covariance function above reads 

2 poo roo /*oo 4 ^^ 4^2 

EX{t)X{t + h) = „ / / e~'^'^^W2~^dw2 / — ;- - ;— 

^ ^ ^ ^ r^P) Jo 1 Vo ' 7-00 22^^ 


1 (j!4-fl)^ 


-dxi 


=4(7^ 


+“ /|xi||a:2-h| 




4^ y 2 • ■' ^"2 

We now pass to the general even-order fractional equation (1.4). 

Theorem 3.2. The representation of a solution to the equation (1.4) can be written as 

roo ^+oo 

X{t) = / w^~^e~^'^ / U2n{x,w)S{t + x)dxdw, P > 0, fi>0. (3.7) 

J. (P) Jo J-QO 

Moreover, the spectral function of (3.7) is 


/(t) = 


a 


[pL + 


and the eovariance function reads 


a 


Covx{h) = -^^[u 2 n{h,W 2 p)] 
where W 213 is a gamma r.v. with parameters /r, 2/3. 

Proof. The solution U 2 n{x,t) to 


d_ _ /_iNn+l d'^'^ 

^ ' (9x2^ 


U2n 


has Fourier transform 
We write 

Since 

we also write 


UiP,t) = = e-y-i. 


e = 




/. 

U{—iP,t) = e 

/-OO g 

/ e^^U2n{x,w) dx 

J —00 


e*^et U2n{x, w) dx. 


In conclusion, we have that 

xit) = (^^ + {-ly 




,2n \ y 


dt‘^^' 


1 r 

W)Jo 


Sit) 


r(/3) 

1 


dwe->^^wJ-^ 


r+oo 

—00 
'+00 


dx U2nix, w) e^StSit) 


m Jo 

and this confirms (3.7). 

From (3.7), in view of (2.20), we obtain 

„2 poo 


EXit)Xit + h) = 


r2(/3) Jo 

+00 


dwi w 


pc 

y-ig-M^'i / 
^0 

/ +00 r+oo 

dxiU2nixi,Wi) / 

-00 J —00 


dw2W2 ^e 


(3.8) 


(3.9) 


(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 


/ H-OO 

dxU 2 nix,w) Sit + x) (3.16) 

-00 


dX2 U2n ix2, W 2 ) 5(3^2 - Xi + h) 
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roo 

/ dw2W2~^e 

Jo 

/ + 0O 

dxi U2n {xi,Wi) U2n {xi - h, W 2 ) 

-00 


^W2 


r2(/3) Jo 

r* + 00 


— 00 

2 roo 


(7 

"'^)J0 

^2 


poo poo 

/ dwi / dw2 tt>2 U2nih, wi + W 2 ) 

Jo Jo 


a 


= ^Eu 2 nih, Wi + W 2 ). 




By following the same arguments as in the previous proof, we get that 

) 

dw U 2 n{h, w) 


EX{t)X{t + h) =^Eu2n{h, W2P) = ^ 


The spectral density of X (t) is therefore 


fix) = 


a 


r(2/3) 




a 




□ 


Theorem 3.2 extends the results of Theorem 3.1 when even-order heat-type equations 
are involved. 

We now pass to the study of the equation (1.9) for n = 1 and n = ^1, 


d^\^ 


(3.17) 


Theorem 3.3. The representation of a solution to the equation (3.17) can he written as 

1 r 
rW)Jo 


^if) ~ wol / ^ ^'^dw ^ f _ Ai ( —) £{t + x) dx dw, /3 > 0, /r > 0. 

J-00 y^w \y 2 >wj 


Moreover, 


and 


fix) = 


a 


{p, -|- KiT^)‘^h 


a 


Covxih) = 


T 


o 




—nh \ 


(3.18) 

(3.19) 

(3.20) 


where Ai{x) is the Airy function and W is a gamma-distributed r.v. with parameters 2/3 
and p. 

Proof. By following the approach adopted above, after some calculation we can write that 


x-{t) = 


r(/3) 


'0 


yjh Ig + T(t) dtc 


is the solution to 


whereas 


dJ\^ 

Xit)=£{t) 


X+(t) = 


1 


r(/3) 


B —1 —LLW — W 

e ^ 


dt^ £{t) dw 


p + ^^J Xit)=£{t) 


(3.21) 

(3.22) 

(3.23) 

(3.24) 


is the solution to 
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The third-order heat type equation 


d 


has solution, for k = — 1, 


u{x,t) = 


:Ai 




X G M, t > 0, 


with Fourier transform 


/ CXD 

-OO 

Formula (3.27) leads to the integral 


(3.25) 

(3.26) 

(3.27) 


e^^u{x,t) dx = e 


te^ 


e E 


because of the asymptotic behaviour of the Airy function (see [1] and [9]). The solution 
to (1.9) with n = 1 (that is k = —1) is therefore (3.21). 

The equation (3.25) has solution, for k = -|-1, given by 


u{x, t) = ^ 7 ^ Ai ( 


—x 


X E M, t > 0. 


Thus, by following the same reasoning as before, we arrive at 


f e^^u{x,t) dx = e 9 

J —OO 


e 


and we obtain that (3.23) solves (3.17) with k = -|-1 is (3.23). 
In light of (2.20) we get 

^2 roo 


E[X-(t) X-{t + h)] 


a 


F2 (/3) 


noc ^ 

J-ooW^ 




-^Swi) -^3102 ' ^ \/3w^ ) 


dwi W 


PC 

r 

Jo 




dw2 w: 


,/3-l 

2 


' —OO 

^2 .OO 


Xi 


;^Ai (dxi 


-n2/m / I e-^^^^dw2wt, 

1 iPj Jo Jo 


PC 

r/ 

Jo 


1 


rAi 




h 


^3{wi + W2) \ ■^3(u;i W2) ^ 




=Ai 


h 


^3^ Vv/3W2^, 

From the Fourier transform (3.27), we get that 


/ ("t) = 


a 


2/3 








=Ai 




dh 


11/2/3 


(/r -|- 

3 

2 —i 2 B arctan — 
(7 e ^ 


(3.28) 


(/l2 + T^)P 
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Also, we obtain that 


E[X+{t)X+{t + h)] 




2/3^ 


=Ai 




with 


/+M = 


a 




2 -\-i23 arctan — 


(/U — {iJ? + r®)^ 


(3.29) 

□ 


Theorem 3.4. The representation of a solution to the following equation 

Q2n+1 \ h 


( 52^+1 \ P 


reads 


roo r-\-oo 

X{t) = / w^~^e~^'^ / M2n+i(NX, /3 > 0, > 0. 

r(p) Jo J-00 

Moreover, the covariance funetion 


a 


Covxih) = -j^Eu2n+l{nh,W2l3)- 

has Fourier transform 

2n + l 

2 2 —i2/3/^ arctan- 

N ^ ^ ^ 

(/i + rar2^+l)2/5 (^2+^2(2n+l))/3 ' 

Proof. The proof follows the same lines as in the previous theorem. 


□ 
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